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• ' The peak time of optical afterglow may be used as a proxy to constrain the Lorentz 

o ' factor r of the gamma-ray burst (GRB) ejecta. We revisit this method by including 

X ■ bursts with optical observations that started when the afterglow flux was already de- 

^ • caying; these bursts can provide useful lower limits on F. Combining all analyzed bursts 

jrt ■ in our sample, we find that the previously reported correlation between F and the burst 

luminosity L^ does not hold. However, the data clearly shows a lower bound Fmin which 
increases with L-y. We suggest an explanation for this feature: explosions with large 
jet luminosities and F < Fmin suffer strong adiabatic cooling before their radiation is 
released at the photosphere; they produce weak bursts, barely detectable with present 
instruments. To test this explanation we examine the effect of adiabatic cooling on 
the GRB location in the L^ — T plane using a Monte Carlo simulation of the GRB 
O ' population. Our results predict detectable on-axis "orphan" afterglows. We also derive 

^_j , upper limits on the density of the ambient medium that decelerates the explosion ejecta. 

We find that the density in many cases is smaller than expected for stellar winds from 
normal Wolf-Rayet progenitors. The burst progenitors may be peculiar massive stars 
with weaker winds or there might exist a mechanism that reduces the stellar wind a few 
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_^_' years before the explosion. 
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1. Introduction 

The prompt emission of gamma-ray bursts (GRBs) is likely produced by dissipative mecha- 
nisms inside the relativistic ejecta of the explosion, while the GRB afterglow is associated with the 



ejecta deceleration by a circum-burst medium (see e.g. |Piranll2004l for a review). The afterglow 
emission is attributed to a relativistic blast wave that involves a pair of shocks - forward and 
reverse. 

One of the most important parameters of GRBs is the Lorentz factor of the relativistic ejecta, 
r, as the models of the prompt and afterglow emissions strongly depend on its val ue. Useful 
constraints on F may be derived using the timescale and spe ctrum of the prompt GRB ( PaczvnskJ 
1986 : iGoodmanlllQSd ). which typically gives F > 100 (e.g. iLithwick fc Saril I2OOII : iHascoet et al.l 
2OI2I ). Another possible way to constrain the GRB Lorentz factor is to use the afterglow peak 
time Tp as a p roxy for the deceleration time of the blas t wave, T^ec, at wh ich the dissipation 
rate peak s (e.g. iMeszaros fc Reed 1 199 71 : 1 Sari fc Piranlll999l ). Using this method iLiang et al.l (|20ld . 
2OI2I ) and lGhirlanda et al.l (120121 ) studied a sample of GRBs with detected optical peaks and found 
a correlation between F and the burst luminosity L^. 

The key assumption of this method, Tp ~ Tdec, is questionable, as optical emission could reach 
its peak at a different time. This possibility is illustrated by the simple model of synchrotron 
emission from a self-similar blast-wave, where the optical light curve can peak at Tp ^ T^cc (e.g. 
Sari et al.lll998l ). However, observations conflict with the late-Tp models and lend some support to 
the Tp ~ Tdcc assumption. In many bursts, optical emission peaks early and steeply, as may be 
expe cted at THprLj The opti c al peak may be dq r ninated by the forward- or reverse-shock emis sion 



(e.g. IMeszaros Sz Reedll997l : ISari fc Piranlll999l : lUhm fc Beloborodov 
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Genet et al.l l2007l ;i . In 



this paper, we accept Tp ~ T^ec as a reasonable assumption and investigate its implications. 

In Sections 2 and 3, we extend the previous analysis by including bursts whose afterglow 
peaked before observations started, which provides a useful upper limit on Tp. We also identify 
the cases where the blast wave at Tp is significantly slower than the ejecta, which corresponds to a 
relativistic reverse shock; in these cases the measurement of Tp provides only a lower bound on F. 
Our analysis does not support the existence of the L^ — F correlation claimed in previous studies. 
Instead, the data shows a lack of bright bursts with low Lorentz factors. In Section 4, we suggest 
an explanation of this fact. Section 5 summarizes our results and discusses implications of observed 
Tp for the nature of the circum-burst medium. 



^This expectation depends on the model of the circum-burst density. A steep rise is firmly pr edicted if the density 
is uniform, p — const, but questionable if p ex R~^ (wind-type medium. IChevalier fc Lj 200C). In the la tter case, 
additional effects such as e pair loading could produce the steep rise toward the peak ( Beloborodovll2002l l. 



2. Sample 
2.1. Bursts with detected optical afterglow peaks 

Table [1] gives the list of GRBs with detected early afte rglow peaks that a r e included in ou r 



sample. This list is a selection from the GRB samples of iLiang et al.1 ([20101) ; iLii et al.1 (J2012l ) ; 



Liang et alj (|2012l ). where we keep only bursts with reliable detection of the peak time. We removed 
(i) GRBs that have optical light-curves with more than one bump, making the peak measurement 
ambiguous, (ii) GRBs for which the peak was measured during a plateau phase (i.e. where the 
optical light-curve is flat in logarithmic scale), and (iii) GRBs with optical light-curves that are 
sampled too sparsely or whose temporal range is too small. 

The redshift-corrected peak times Tp = T° /{I + z) are shown in Figure [1] versus the GRB 
isotropic equivalent gamma-ray energy E^ and the corresponding average luminosity 

L, = ^, (1) 

with 

T - -^ (2) 

where Tgo is an approximate measure of the observed burst duration (time during which 90% of 
the emission is received). One can notice that there is a good correlation between Tp and L^ (or 
E^); we argue below that this correlation is spurious. 



2.2. Bursts w^ith upper limits on Tp 

For some bursts the peak is not observed because observations start too late. These bursts 
are also useful for our purposes, as some of them give strong upper- limits on Tp. A sample of such 
GRBs is listed in Table [2] and the corresponding upper limits are shown in Figure [H The number 
of available strong limits (24) is comparable to the number of peak detections (20). The limits are 
robust. For many bursts in the sample, the optical decay was already well established when the 
observations started, without any evidence for an increasing decay index. This suggests that the 
peak was reached well before the beginning of observations. 

The obtained limits on Tp are never below a few tens of seconds, which reflects the typical 
delay in response of robotic optical telescopes to alerts from 7-ray telescopeso Note also that the 
limits tend to be less constraining for weak (low E^) bursts. This tendency can be explained by 
selection effects. The localization of weak bursts by 7-ray telescopes is more difficult. Weak bursts 
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Due to a fortun ate chain of events, optical observations of GRB 080319B started before the 7-ray trigger 



( Racusin et al.l 120081 ) . But in this special case the rise of the optical afterglow is hidden by the bright prompt 



optical component. 



Table 1: Bursts with detected optical peaks. 



GRB 



E. 



7,52 



T-iobs 
P 



^90 



References 



990123 


1.60 


436.52±60.31 


47±10 


63.3±0.3 


1 


2 


050820A 


2.612 


159.2±12.4 


477±6 


600±50 


2 


3 


060418 


1.489 


48.6±10.6 


170±5 


52±1 


2 


3 


060605 


3.78 


2.8±0.5 


590±45 


19±1 


2 


3 


060607A 


3.082 


23.4±1.5 


179±3 


100±5 


2 


3 


061007 


1.261 


421±41.9 


77±1 


75±5 


2 


3 


070318 


0.836 


1.3±0.3 


507±46 


63±5 


2 


3 


070419A 


0.97 


0.2±0.02 


765±30 


112±2 


2 


3 


071010B 


0.947 


1.7±0.9 


287±145 


35.74±0.5 


2 


3 


071031 


2.692 


3.9±0.6 


1213±2 


180±10 


4 


3 


080603A 


1.67842 


1.5±0.5 


1044±167 


150±10 


5 


3 


080710 


0.845 


0.8±0.4 


1934±46 


120±17 


2 


3 


080810 


3.35 


30±20 


117±2 


108±5 


2 


3 


081008 


1.967 


2.8±0.5 


163±2 


185±39 


6 


3 


08 1203 A 


2.1 


17±4 


295±2 


294±71 


7 


3 


090313 


3.375 


4.6±0.5 


1315±109 


78±19 


2 


3 


090812 


2.452 


45.9±6 


71±8 


70±5 


8 


3 


091029 


2.752 


7.4±074 


328±50 


39.2±5 


9 


2 


100906A 


1.727 


33.4±3 


101±4 


114.4±1.6 


2 


3 


110205A 


2.22 


56±6 


948±3 


257±25 


2 


3 
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also have fainter afterglows (JGehrels et al.ll2008l ). making their early observation/detection more 



challenging. The tendency of missing the peak for weak GRBs is enhanced by the fact that they 
are preferentially detected at smaller redshifts; then the time dilation factor (1 -\- z) is less efficient 
in moving the peak to a later observational time. While it is easy to miss an early optical peak, 
Tp < 100 s, we are not aware of any selection effects that could lead to preferential non-detection 
of late peaks Tp ~ 10^ — 10^ s. In this range, the data should represent the true distribution of Tp. 

The data presented in Figured] may be summarized as follows: there is no intrinsic correlation 
between L^ (or E^) and Tp. Instead we observe a lack of bright bursts with late afterglow peaks. 
For a given L^, there appears to exist a maximum peak time Tp^maxlT-y) which corresponds to the 
blue boundary in Figure [TJ A crude approximation to this boundary is given by Tp^max(-^7) '^ 
200 L-%' s. 



3. Estimates for F 

The GRB afterglow is likely emitted by the blast wave resulting from the interaction of the 
relativistic ejecta with the ambient medium. The blast wave involves two shocks: the forward shock 
sweeping the external medium and the reverse shock propagating back into the relativistic ejecta. 
As discussed in Section [H it is reasonable to assume that the afterglow peaks at the deceleration 
time T(jec when most of the ejecta energy has been transmitted to the blast-wave through the 
reverse shock. This happens at the "deceleration radius," 

-^dec = -r— Tt.2 !■" ' (3) 



V 4vrc2 F^^pde 

where Ef.^ is the energy of the ejecta and Fbw is the Lorentz factor of the blast wave at Rdec] 
Pdec is the external density at -Rdco and s describes the slope of the external density profile, 
s = dlnp/dlni?. A uniform medium is described by s = and a wind medium by s = 2. We 
assume that the afterglow peaks at the deceleration time, 

TpC^Tdec--^. (4) 

Then the measured Tp provides an estimate for the blast wave Lorentz factor, 

1 
3 — s 1 — r] E^ 



bw 



327rc5 rj pdecT^ 



(5) 



where r] is the fraction of the initial energy of the GRB ejecta that is converted into prompt 
radiation. 

Two aspects of Equation ^ should be noted: (1) the estimate depends on the poorly known 
ambient density as [/Jdec/(3 — s)]^ and the prompt efficiency as [ri/{r] — 1)]^ ' . (2) The estimate 
gives the Lorentz factor of the blast wave, Fbwj not the ejecta Lorentz factor F. 



We expect F ~ Fbw if Tp ^ T-y. Indeed, Tp is associated with the time it takes the ejecta 
to transfer most of its energy to the blast wave, i.e. the time it takes the reverse shock to cross 
the main, most energetic part of the ejecta of thickness A ~ cT^. The crossing time A/v^s is long 
and gives Tp ~ {c/vrs)T^ 3> T^ if the reverse shock is non relativistic, v^s "C c, which is equivalent 
to r ~ Tbw In this case Equation ([5]) effectively gives an estimate of the mean Lorentz factor F 
of the ejecta. In contrast if Tp < T^, the reverse shock may be highly relativistic. Then Fbw is 
significantly smaller than F and Equation ([5]) significantly underestimates F. GRBs with Tp < T^ 
are highlighted in red in Figure [TJ 

Note that T^ that we use as a measure of the GRB duration may overestimate the duration of 
the main part of the GRB if the burst has a temporally extended tail of relatively weak emission. 
For such bursts Tdcc < T^ is possible. A better estimate for T-y would give T^ec — T^, so that Tp is 
not smaller than T-y. In agreement with theoretical expectations, we found no burst in our sample 
where the optical afterglow peaks before the main part of the GRB emission has been received. 

Equation ([5]) assumes a static external medium and neglects the fact t hat the prompt GRB 



radiation exerts pressure and accelerates the medium ahead of the blast wave (JThompson &: Madau 



2000l : lBeloborodovll2002l ). This pre-acceleration is strong (relativistic) up to the radius, 

R,cc = 2xlO^''iE^,52)^/^ cm. (6) 

If Race exceeds Rdec given by Equation (J5]), the true deceleration radius is increased, and the 
dissipation rate peaks at 

^p = |Fj = 4(^.-)"^(^)".. (7) 

We checked for our GRB sample that i?acc < ^dcc as long Pdcc/'^p ^ 10^ cm~^. For most bursts 
in our sample, /9dec/"^p ^ 10^ cm~^ would imply low values for F, which would contradict the 
constraint from the prompt emission (see discussion in Section 15. 2p . Therefore we only consider 
Pdec/"ip ^ 10'^ cm^^ and neglect the pre-acceleration effect for our sample. 

Using Equation ([5]) we estimated F for each burst in the sample (Figure [2]). We fixed r] = 0.5 
and s = 2 in our numerical estimates; the uncertainty in their exact values weakly affects the 
results. If only the peak "detections" are considered. Figure [2] would suggest that there is a 
correlation between F and L^. However, the numerous lower limits show that the upper part of the 
diagram must be broadly populated by GRBs. Also note that four detections of Tp correspond to 
the relativistic reverse shock regime and give only lower limits on F. We conclude that there is no 
evidence for a correlation between L^ and F. However the lack of bright bursts with low F appears 
to be a robust feature. 

Figure m illustrates how the GRB sample is transformed from L^ — Tp plane to L^ — Fbw plane, 
following the relation between Fbw and Tp, Fbw oc E-y Tp (Equation ([5])). This transformation 
compresses the sample along the Fbw axis and induces a dependence of Fbw on L^ with a positive 
slope of 1/8. In combination with the selection effect that suppresses short-Tp and low-L^ bursts 
in the sample, this enhances the spurious correlation between L^ and F. 



As mentioned in Section [H the key assumption that the optical afterglow peaks at the decel- 
eration radius may not be reliable. Therefore, it is useful to consider a more general model where 
the afterglow peaks at a radius Rp related to F and L-^ by 

Rp oc r"L^ . (8) 

The simplest model with p^ec = const corresponds to a = —2/3 and /3 = 1/3 (see Equation Q). 
For other values of a, /?, the observed upper boundary T™^^ oc L^ (where A ~ —3/5; see Figure [I]) 
in the L^ — Tp diagram still transforms into a lower-boundary Fmin oc L-,, in the L^ — T 

diagram. As an illustration, Figure [3] shows F estimated assuming a = /3 = (i.e. Rp = const). 
The results are similar to those in Figure [5J The conclusion that Fmin grows with L^ holds as long 
as (A — 0)/{o: — 2) > 0. Violation of this condition would require a and /3 that are significantly 
differential from a = —2/3 and /3 = 1/3. If a ~ 2 then Tp weakly depends on F, so it cannot be 
used to estimate F. 



4. Lack of bright bursts with a large Tp-. a result of adiabatic cooling? 

4.1. Adiabatic cooling below the photosphere 

The photospheric radius i?* (the radius where the explosion ejecta becomes transparent to 
Thomson scattering) is given by 

R^ ~ 2.9 X 10i2^o.2i52F^^ cm , (9) 

where L is the isotropic power of the flow, and k is the Thomson opacity (in units of 0.2 cm^g^^); 
K may be significantly increased by pair creation. For the most powerful explosions (L^ ~ 10^^ — 
-j^q54 gj,g g-i^ with relatively low Lorentz factors (F < 10^) the photospheric radius is exceptionally 
large. If the GRB radiation is produced at a smaller radius -Rdiss (where internal dissipation peaks) 
the burst may be buried by the large optical depth, since it implies strong adiabatic cooling of the 
radiation trapped in the expanding ejecta. 

To illustrate this possibility, suppose that the GRB emission is generated before the ejecta 
reaches a radius 

-Rdiss = r -Ro 5 (10) 

where Rq is a fixed constant. The scaling of the cutoff radius -Rdiss oc F^ is expected for mechanisms 
that dissipate the energy of internal motions or magnetic energy. Combining Equations ([9]) and 
(fTO|) . one finds the optical depth at i?diss) 

Tdiss = iP^ OC LF-5 . (11) 

-Kdiss 



■^This would be possible e.g. if bright bursts have much denser ambient medium with dlnpdcc/dlnL-y > 3. 



If the dissipation occurs far below the photosphere (rdiss ^ 1), the resulting radiat ion released 



-2/3 



Then for 



at the photosphere is adiabatically cooled by the factor of 2t^^J (jBeloborodovl 12011 

given ejecta power L and dissipation efficiency at -Rdiss the observed burst luminosity scales as 



L^ oc rio/3 . 



(12) 



The slope of this relation is suggestively close to the slope of Fmin — L-y relation seen in Figure [2j 



4.2. Properties of adiabatically cooled bursts 

Even though the dissipation mechanism at r < i?diss can be non-thermal, the produced radia- 
tion will be progressively thermalized during the subsequent adiabatic expansion between i?diss 
and i?*. Let T^ be the electron temperature (measured in the ejecta frame), and Tdiss > 1 
be the Thomson optical depth at i?diss- Two different regimes can be distinguished: (i) For 
■^diss ^ meC^/SkTf, ~ 511 (3A;Te/l keV)~ , an exponential cutoff must form in the radiation spec- 
trum at energy 



TnieC^ 



F 



^max ~ — ~ 511 keV . (13) 

''"diss ''"diss 

This cutoff is a result of significant Compton downscattering (recoil effect) at E > -Bmax! the 
spectrum at £' < -Emax is weakly affected, (ii) For Tdiss ^ rrieC^ /"ikT^, the spectrum is exponentially 
suppressed above SfeTeF. The low-energy part of the spectrum is also affected by multiple Compton 
scattering — the spectral slope steepens as the photons tend to thermalize with electrons. 

Thus, one expects significant chang es in the burst s pectr um a fter strong a diabatic cooling. 
Such unusual GRBs have been observed. iGhirlanda et al.l (J2003l ) and iRvdd (J2004l ) found that some 
bursts have spectra with very hard low energy indices and possibly exponential cutoffs. Similar 
quasi-thermal GRBs are found in both BATSE (Kaneko et al. 2006) and Fermi GBM (Goldstein 
et al. 2012) catalogs. These bursts — especially those with a low peak energy — may be generated 
by adiabatically cooled explosions with a large photospheric radius i?*. 

Another expected feature of bursts with small F and large i?* is the suppression of variability 
on short timescales. In these bursts, the minimum variability timescale Atobs ~ i?*/2F^c can be as 
large as 10 s (see Equation pT|) below), and their lightcurves are expected to be smooth. 

While adiabatic cooling can significantly reduce the emitted GRB energy E^, the ejecta energy 
E'ej remains large. Thus, adiabatically cooled bursts are expected to have unproportionally bright 
afterglows. They should lie in the upper part of the gamma-ray fluence/X-ray afterglow flux 
distribution, which spreads over two orders of magnitude (see for example Fig. 2 in iGehrels et al.l 
20081 ). As we show below, the majority of GRBs with strong adiab atic cooling avoid d etection. 



Then they become prime candidates for "orphan" afterglows (e.g. iHuang et al.l |2002| ). These 



afterglows are expected to peak at late times Tp, as they are generated by the low-F ejecta. 



4.3. Monte Carlo simulation of a GRB population 

To illustrate how bright bursts with low Lorentz factors are depleted by adiabatic cooling, we 
produced a synthetic GRB population using a Monte Carlo simulation, with the following assump- 
tions: 



(1) We adopted the GRB rate Rsfs from lPorciani fc Madaul (12001). We assume that the rate keeps 



increasing at z > 2 as suggested by observations (e.g. iDaigne et al. 
Salvaterra et al.ll2012l ): we cut it off at Zmax — 20. 



200d : IWanderman fc Piranlboiol : 



(2) The GR B luminosity functio n is assumed to follow a power-law distribution of index —1.5 as 
suggested by lOaigne et al.l (J2006l ). in a broad range of 10^° erg s^^ < Lj < 10^^ erg s~^. 

(3) For each GRB, the spectrum of generated radiation at i?diss is assumed to be a broken power-law 
with a low energy index a = — 1 and a high energy index /3 = —2.5. The peak energy of the spec- 
trum, Ep, is assumed to correl ate with L^; we u se the relation Ep ~ 300(L^/10^^ erg s~^)^/^ keV 
with a scatter Cdex = 0.3 (e.g. iNava et al.ll2012l ). 

(4) The radiation spectrum is injected with r] = 0.5 at i?diss = 6 x 10^^ (F/lOO)^ cm (which would 
correspond to a variability timescale Aiobs = ^diss/2r^c ~ 10 ms if i?diss > R*)- 

(5) The logarithm of Lorentz factor T is randomly chosen for each burst from a uniform distribution 
in the range 1 < logF < 3. 



The adiabatic cooling effect is calculated as follows: 

(6) The photospheric radius of each GRB is obtained from Equation Q, with ko.2 = 1- If 
-Rdiss < 2~'^/^i?*, the burst is cooled by a factor of 2(i?*/i?diss)~^ changing L^ and Ep from 
their initial values at .Rdiss- The burst spectrum is changed as explained in Section 142} the cutoff 
at -Emax is approximated by a step function. 

(7) The simulated GRB is assumed to be detected if its observed photon flux in the Swift band 
15 - 150 keV is above a threshold of 0.2 ph cm'^ s"^ (JBandlbood ). 



The results of our simulation are shown in Figure [5j One can see that adiabatic cooling depletes 
the low-F/high-L^ corner of the L^ — F diagram. The resulting distribution resembles the observed 
one in Figure [2j This simulation also allows one to estimate the impact of adiabatic cooling on 
the observed distribution of Ep (Figure E]) . A burst that suffers adiabatic cooling is moved along 
a track Ep oc L^, which tends to create GRBs with Ep below the original correlation. However, 
the effect on the population of detected GRBs is weak, because cooled bursts become undetected 



(due to the reduced Ep and the spectral cutoff at E^ 
correlation. 



before they become outliers in the Ep — L^ 
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Discussion 



5.1. Adiabatic cooling 

If the peak time of optical afterglow, Tp, is indeed a good proxy for the blast wave deceleration 
radius -Rdec, observations imply a lack of bright bursts with low Lorentz factors (Section [3|). We 
argued that this lack may be expected, as energetic explosions with low Lorentz factors should have 
unusually large photospheres iZ* > 10^^ cm. Since the dissipation mechanism generating radiation 
in these bursts is likely limited to smaller radii, i?diss ^ R*, the burst is expected to suffer strong 
adiabatic cooling and become undetectable (Section 4). 

Most models of the prompt GRB emission place the emission source at radii smaller than 
10^^ cm, especially if F is small. For example, dissipation of internal motions or magnetic energy 
in the ejecta is expected to end at a radius that scales as F^ and becomes smaller than R^ at 
small F. Then the trapped radiation is adiabatically cooled and the burst becomes inefficient. This 
argument is applicable t o any dissipation mechanism ge nerating the burst — e.g. collisionless sh ocks 
(JRees Sz Meszaroslll994: iDaigne Sz Mochkovitchlll998l ). collisional heating (JBeloborodovl l20ld ) . or 



magnetic reconnection (jSpruit et al. 



2001 



). Note also that the neutron component of the jet, which 
can play a significant role in collisional dissipation, may not survive to i?=K in low-F bursts. The 
mean radius of neutron decay is Rp ~ 9 x 10^^(F/30) cm, and its ratio to the photospheric radius 
is given by 



:^ ^ 7.5 X 10- 



«:o.2 L 



54 



30 I «'■ 



(14) 



The cooled bursts still produce energetic ejecta that can drive an energetic blast wave in 
the external medium and generate bright afterglow emission. We argued that one could observe 
"orphan" afterglows from such explosions, with undetected prompt GRBs, even when the burst 
is observed "on-axis," i.e. the relativistic jet is directed toward the observer. As this paper was 
completed, Paloma r Transient Factory detected an event consistent with on-axis orphan afterglow 
(|Cenko et ahlboid ). 



We also argued that the prompt emission of strongly cooled bursts can be occasionally detected. 
As discussed in Section 4, these bursts have special properties. They should have soft spectra 
resembling quasi-thermal emission and their light curves should be smooth. 



5.2. GRB ambient medium: a low density wind? 



As some long GRBs are associated with Type Ib,c supernovae, their progenitors are expected 
to be Wolf-Rayet stars. Wolf-Rayet stars in our galaxy are observed to los e mass at a typi cal rate 
M ~ 1O~^M0 yr~^ through strong winds of velocity w ~ 10® cm s~^ (e.g. ICrowtherl 120071 ) . If the 
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wind has constant M, its density follows the R ^ density profile, 



p{R) 



where 



A 



M 



M 



4ttw \ W-^Mq yr-i 



Ws 



A 



-' Ao, 



^0 = 5 X 10^^ g cm^^ . 



(15) 
(16) 



As discussed in Section [3l the peak time Tp of the optical afterglow provides an estimate of 



the ejecta Lorentz factor T ~ Fbw if Tp ':$> T^. This estimate scales as p^^^ . On the other hand, 
an independent constraint on T can be derived by considering the Thomson opacity of the ejecta 
and its photosphereo The photospheric radius R^: (Equation Q) implies a minimum variability 
timescale Atobs ~ i?*/2r^c. Thus, the observed Atobs is expected to satisfy the condition, 



At, 



obs 



> 



R* 



1 + z ^ 2r2c 



5^0.2- 



'-L 



V 



7,52i 2 



ms . 



(17) 



From this condition, using the observed L^ and Atobs > one obtains a lower bound on the ejecta 
Lorentz factor. Then, combining Equations ([5|) and (fT7|) . an upper limit on p^ec can be derived. 
For a wind medium this limit translates into an upper bound on A, 



A<A, 



10"ko.2-^/^ 



rj 



rj 



1/5 



P 1/5 
-C/7,52 



At, 



obs 



l + Z 



4/5 



r//5Tp,2-i gcm-i 



(18) 



This upper limit is sensitive to the Thomson opacity k, which may be significantly increased by 
e creation in the ejecta. The numerical value for ^max in Equation (jlSp is given for the lowest 
possible K = 0.2 cm^ g~^, which gives a conservative upper limit on A. Figure [7] shows the upper 
limits obtained for our GRB sample, where we used rj = 0.5, K0.2 = !> and Atobs = 1 s. In some 
bursts, the inferred ^max is well below Aq. For GRB 060605, ^max is ~ 6 times smaller than ^q- 
Even stronger limits could be derived with more sensitive detectors and more detailed analysis of 
variability timescales. GRBs that give the lowest Amax have relatively low luminosities, and poor 
photon statistics make it difficult to see the true minimum Atobs, which might be shorter than 1 s. 

Figure [7] suggests that at least some GRB progenitors are peculiar massive stars whose winds 
are weaker than typical Wolf-R a.yet stars observ ed in our galaxy. This may be the result of a 
lower metallicity of the star (e.g. IVink et al.ll200ll ). a property that se ems to be preferre d by GRB 



progenitors, and agrees with observations of their host galaxies (e.g. iPerley et al.ll2013l ). Weaker 
winds extract less angular momentum from the progenitor, leading to coll apse with faster rotation , 
which is required for GRB central engines producing collimated jets (e.g. IWooslev &: Hegeiil2006l ). 



Note also that GRBs with bright optical afterglows (such as bursts in our sample) imply a selection 



'For GRBs with detected high-energy emission (above 100 MeV) another co nstraint on F could be der ived from 



the requirement th at the high-energy photons avoid 7-7 absorption (e.g. iLithwick fc Sari 
Hascoet et al.ll2012h . There are no such bursts in our sample. 



2001 



Granot et al.ii2008i : 
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against high-metallicity environment, which tends to obscure the optical emission (e.g. iLevesque 



20131 ). 



The low wind density could also be explained by a change in the stellar mass loss rate M 
shortly before the explosion. The wind medium at a characteristic radius i?dec ~ 10^^ cm was 
ejected by the progenitor R^ec/w ~ 3 yr before the explosion. The mass loss rates of Wolf-Rayet 
stars in the last few years of their lives are uncertain; M might decrease as the star evolves toward 
the collapse. 

This work was supported by NSF grant AST- 1008334. 
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Table 2: Bursts with upper limits on Tp 
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Fig. 1. — Bursts on Tp — E^ and Tp — L^ planes. Crosses represent bursts with detected Tp, and 
arrows show upper hniits. Bursts with Tp < T^ are shown by red crosses and arrows. The blue line 
shows the observed boundary of the burst population. 
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Fig. 2. — Estimated Lorentz factor of the GRB blast wave, Tbw) at the deceleration radius. The 
product Tbw X Pdec^ is estimated using Equation ([5]) (with rj = 0.5) and shown versus the burst 
luminosity L^ (Equation ^). The ambient density at the deceleration radius, pdec, is normalized 
to po/rup = 1 cm^^. Bursts with Tp < T^ are highlighted in red color; for these bursts the ejecta 
Lorentz factor F can be substantially higher than Fbw (see text). 
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Fig. 3. — Estimated Lorentz factor of the GRB blast wave Fbw = (-Rp/2cTp) ' for an arbitrarily 

— 1/2 

fixed afterglow peak radius Rp. The product Fbw^p shown in this figure is proportional to 



-,-1/2 



SO the diagram is a simple transformation of the Tp — L^ diagram shown in Figure [TJ 
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Fig. 4. — Effects of the Tp — )• F transformation. This schematic figure illustrates how the GRB 
population presented in Figure [U (Tp — L^ plane) transforms into the V — L^ plane. The hatched 
and colored areas are indicated to better visualize how different regions on the Tp — L^ plane 
transform to the T — L-y plane. The red (horizontally hatched) region is where strong selection 
effects are expected to suppress the observed population. Combined with the real lack of bursts in 
the blue (vertically hatched) region, this leads to a spurious correlation between Tp and T^. The 
corresponding spurious F — L^ correlation is enhanced by the deformation of the population in the 



F — L-y coordinates 
in the lower panel. 



the black rectangular in the upper panel is transformed into a "parallelogram" 
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Fig. 5. — Effect of adiabatic cooling on a GRB population in the L^ — T plane. Left: observed 
distribution of bursts when adiabatic cooling is not included. Right: observed distribution of bursts 
when adiabatic cooling is included. Red crosses represent bursts that have suffered cooling and 
remained detectable. Many more cooled bursts became undetectable and disappeared from the 
diagram. The circled red crosses show detected bursts with large cooling factors / > 10 (1.3% of 
detected bursts). These bursts are expected to have special spectra of the prompt GRB emission 
(see text). 
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Fig. 6. — Effect of adiabatic cooling on a GRB population in the L^ — Ep plane. The simulation 
is the same as in Figure [5l Left: observed distribution of bursts when adiabatic cooling is not 
included. Right: observed distribution of bursts when adiabatic cooling is included. Red crosses 
represent bursts that have suffered cooling. The circled red crosses show detected bursts with large 
cooling factors / > 10 (1.3% of detected bursts). Strong adiabatic cooling can in principle create 
outliers of the initial Ep — L^ correlation, however most of them become undetectable and disappear 
from the diagram. 
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Fig. 7. — Constraints on the wind from the GRB progenitor. The figure shows the derived upper 
hmits on the wind density parameter A = M/Attw in units of ^o = 5 x 10^"^ g cm~^. 



